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ABSTRACT 

The X-ray spectra of Galactic binary systems dominated by a quasi-thermal com- 
ponent (disc dominated or high/soft state) are well described by a standard Shakura- 
Sunyaev disc structure down to the last stable orbit around the black hole. This is not 
the case in the very high (or steep power law) state, where the X-ray spectra show both 
a strong disc component and strong, steep tail to higher energies. We use simultaneous 
optical- ASCA-RXTE data from the black hole transient XTE J1550 — 564 as a specific 
example of two such spectra, where the power emitted in the tail is more than 50 per 
cent of the bolometric luminosity. These have disc spectra which are significantly lower 
in temperature than those seen from the same source at the same luminosity in the 
high/soft state. If these give a true picture of the disc then either the disc emissivity 
has reduced, and/or the disc truncates above the last stable orbit. 

However, it is often assumed that the tail is produced by Compton scattering, 
in which case its shape in these spectra requires that the Comptonising region is 
marginally optically thick (r ^ 2-3), and covers a large fraction of the inner disc. This 
will distort our view of the disc, especially of the hottest temperature material. We 
build a theoretical model of a Comptonising corona over an inner disc, and fit this to 
the data, but find that it still requires a large increase in inner disc radius for a standard 
disc emissivity. Instead it seems more probable that the disc emissivity changes in 
the presence of the corona. We implement the specific inner disc-corona coupling 
model of Svensson & Zdziarski (1994), in which some fraction / of the accretion 
power is dissipated in the corona leaving only a fraction 1 — / to be dissipated in 
the optically thick disc. We show that this can explain the low temperature/high 
luminosity disc emission seen in the very high state with only a small increase in radius 
of the disc. While this inferred disc truncation is probably not significant given the 
model uncertainties, it is consistent with the low frequency QPO and gives continuity 
of properties with the low/hard state spectra. 
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1 INTRODUCTION 

Galactic black hole binaries at high mass accretion rates 
show spectra which are often dominated by a soft, quasi- 
thermal component (e.g. the review by Tanaka & Lewin 
1995). This shape matches very well with that expected 
from an optically thick, geometrically thin accretion disc, 
as described in the seminal paper by Shakura & Sunyaev 
(1973). This identification is strengthened by observations 



of variability in which the temperature, Ti„, and luminosity, 
L, of this quasi-thermal component change together in such 
a way as to indicate that the size of the emitting structure 
remains approximately constant (L oc T^„)- The last stable 
orbit is the only obvious candidate for such a radius, so this 
gives compelling evidence for a key prediction of general rel- 
ativity, as well as providing an observable diagnostic of the 
mass and spin of the black hole (Ebisawa et al 1994; Kubota 
et al 2001; Gierliiiski & Done 2004; Kubota & Done 2004; 
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hereafter KD04). The most recent models of accretion disc 
spectra are able to reproduce this behavior (Davis et al 2005; 
Davis, Done & Blaes 2006). These supersede earlier models 
which showed strong deviations from the L oc T*^ relation, 
but which did not include metal opacity nor treat the dif- 
ference in vertical structure between the radiation and gas 
pressure dominated regimes (Merloni, Fabian & Ross 2000). 

The disc dominated spectra always contain some small 
fraction of emission in a tail extending to much higher en- 
ergies. However, some high luminosity spectra have a rather 
different spectral shape, where the tail contains a much 
larger fraction of the bolometric luminosity. The disc emis- 
sion in these spectra (variously termed very high or steep 
power law dominant state, hereafter VHS; Miyamoto ct al. 
1991; McClintock & Remillard 2003) is less obvious due to 
the lower contrast against the strong continuum. Nonethe- 
less, it can be identified in spectral fitting, but its derived 
parameters (temperature and luminosity) are somewhat de- 
pendent on how the tail is modelled, which in turn depends 
on the radiation processes involved. A power law or bro- 
ken power law description for the tail might be appropriate 
if the emission is from nonthermal electrons in the jet, but 
these models require fine tuning in order to match all the ob- 
servational constraints. Conversely, Comptonisation models, 
where the disc is the source of seed photons, rather natu- 
rally reproduce both multiwavelength spectral constraints 
and X-ray variability patterns (see §3), so here we explore 
the consequences of such models. 

Firstly, the disc temperature in GBH is close to (or 
within!) the observed X-ray bandpass, so the low energy 
break in the Comptonised component is predicted to be 
present, giving a rather different spectrum at low energies 
to that of a power law which extends unbroken through the 
disc emission. Thus using proper Comptonisation models 
can significantly change the parameters of the inferred disc 
component (Done, Zycki & Smith 2002) and distort the mea- 
sured luminosity-temperature relation (Kubota ct al 2001; 
Kubota & Makishima 2004; KD04). This helps considerably 
in interpreting data from the weak VHS, where the disc con- 
tributes more than half of the total bolometric luminosity. 
A power law model for the tail in these data gives derived 
disc temperatures which are significantly higher than those 
seen from the same source at the same disc luminosity in 
the high/soft state. Using proper Comptonisation models in- 
stead of the power law approximately recovers the L oc Ti„ 
behavior (Kubota ct al 2001; Kubota & Makishima 2004; 
KD04). The match to the high/soft state data is even better 
with an additional correction to the observed disc luminos- 
ity to account for the photons which are scattered into the 
Comptonised spectrum (Kubota & Makishima 2004; KD04). 

However, the most strongly Comptonised VHS (where 
the disc is less than half of the bolometric luminosity) have 
disk spectra which show significantly lower temperatures 
than those seen at the same disc luminosity in the high/soft 
state. These are inconsistent with a constant disc radius even 
with proper Comptonisation models and correction for scat- 
tering (KD04). This could indicate that the disc is trun- 
cated above the last stable orbit (KD04), but here we ex- 
plore whether this conclusion is required with a more sophis- 
ticated approach to Comptonisation. The parameters de- 
rived for the hot electrons indicate that the Compton cloud 
is marginally optically thick, with r ~ 2-3. This material 



plainly must intercept the disc photons in order to Comp- 
tonise them, so will also distort our view of the disc (KD04). 
An optically thick corona over the inner disc could simply 
hide the hottest part of the disc from view. More likely, the 
power dissipated in the corona can reduce the disc emissiv- 
ity (e.g. Svensson & Zdziarski 1994) since the energy in the 
corona must derive ultimately from the accretion flow. We 
build explicit models of an inner disc- corona to explore these 
ideas, and fit them to two strong VHS spectra (simultane- 
ous 0.7- 200 kcV ASCA and RXTE data) of the Galactic 
black hole binary XTE J 1550 — 564 to constrain the possible 
geometries of the disk-corona in this state. 

In §2, we review the spectral parameters (luminosity 
and temperature) of the geometrically thin disk in the disk 
dominated high/soft state. In §3 we describe the simulta- 
neous ASCA-RXTE observations of the broadband spectra. 
We compare the derived disc emission in the VHS to that 
seen in the high/soft state, showing that it implies a very 
different disc area and/or emissivity (§4). In §5 we describe 
our model of an inner disc corona above a standard emissiv- 
ity disc. Fitting this to the data shows explicitly that optical 
depth effects alone cannot hide the inner disc, and that this 
model requires a large change in disc area compared to the 
high/soft state. We modify this in §6 to incorporate ener- 
getic coupling of the inner disc and corona as in Svensson 
& Zdziarski (1994). The reduction in disc emissivity and 
temperature under the corona means this model can fit the 
observed spectra with only a small - or possibly no - increase 
in disc area. 

We review all the results of the models in §7, and then 
incorporate additional constraints from other spectral and 
timing characteristics of this state in §8. We show that an 
optically thick Comptonisation in an inner disc corona is 
not inconsistent with the presence of high frequency QPO's. 
However, it does conflict with claims of a strong, extremely 
broad iron line in these spectra (Miller ct al 2004) . Wc show 
that there is strong evidence supporting Comptonisation as 
the emission mechanism for the tail (§9), whereas there is an 
alternative explanation for the iron line, in which its width 
is overestimated due to the unrecognized presence of absorp- 
tion lines (Done & Gierlinski 2006). 

All the Comptonisation models considered here require 
that the inner disc changes in extent and/or emissivity in 
response to the presence of the corona. There are no solu- 
tions in which the inner disc has the same radial structure 
as that seen in the high/soft state. 



2 THE DISK DOMINANT HIGH/SOFT STATE 
AS A TEMPLATE FOR THE VHS 
GEOMETRY 

In all the following we use the binary system parameters 
determined by the optical studies of Orosz et al. (2002) i.e. 
mass of 8.4-11.2 Mq, distance, D, of ~ 5 kpc and binary 
inclination angle, i, of 70° (Orosz et al. 2002). However, 
we are primarily interested in changes in the disc emitting 
area rather than absolute values i.e. in comparing the disc 
structure in the VHS with that seen in the disc dominated 
high/soft states. The high/soft states in XTE J 1550 - 564 
have been well fit using the diskbb model (Mitsuda et al. 
1984) in XSPEC (Kubota & Makishima 2004; Gierlinksi & 
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Done 2004; KD04), giving a clear L oc T*'^ relationship, im- 
plying constant area, over a factor 100 change in disc lumi- 
nosity. Such an obviously thermal relationship between lu- 
minosity and temperature from DISKBB fits is now predicted 
by the best current models of the disc spectra which include 
radiative transfer with full metal opacities and self consis- 
tent vertical structure as well as relativistic effects (Davies 
et al 2005). 

The apparent inner disc radius of this source, rin, is ob- 
tained from the normalization of diskbb model, defined as 
r:„ cos i/(D/10kpc)^ The average value inferred from all the 
high/soft state constant area diskbb results is 59 km. How- 
ever, the 'true' inner disc radius is related to this apparent 
radius via a series of correction factors (for the stress free in- 
ner boundary condition, effective temperature not equal to 
colour temperature and relativistic effects; e.g. Kubota et al 
1998; 2001). However, the fractional change in inferred disc 
radius is unaffected by uncertainties in distance, inclination 
and correction factors, so we give the apparent inner radius, 
rin as a ratio of diskbb normalisation seen in these fits to 
that derived from all the high/soft state data in Kubota 
& Makishima (2004). These ratios are not sensitive to any 
of the binary system uncertainties, so can sensitively show 
the change in the disc properties between the high/soft and 
VHS. 



3 ASCA-RXTE SIMULTANEOUS 
OBSERVATIONS 

ASCA observations of this source were performed three 
times, on 1998 September 12, 23, and 1999 March 17. There 
are simultaneous RXTE observations corresponding to each 
ASCA observation. The second data set shows a strongly 
Comptonised VHS spectrum, and was included in the anal- 
ysis of KD04. Here we also study the first simultaneous 
data set, which is at the boundary between the low/hard 
state and the strongly Comptonised VHS. We do not use 
the third ASCA dataset as the gain was changed on the 
corresponding RXTE data dataset midway through the ob- 
servation, and the spectrum is only a weakly comptonised 
VHS. The observational ID's of the corresponding RXTE 
data are 30188-06-05-00 and 30191-01-10-00, for the first and 
second ASCA observations, respectively. We use the data as 
extracted by KD04 (standard procedures for selecting good 
PCA and HEXTE data, with PCA responses derived for 
each observation using pcarsp v8.0: see Appendix for a com- 
parison with pcarsp vlO.l). We use the PCA data from 3- 
20 keV, and HEXTE from 20-200 keV. 

The ASCA CIS events were extracted from circular re- 
gions of 6' radius centered on the image peak, after selecting 
good time intervals in a standard procedure. Net exposures 
of 1.4 ks and 2.4 ks were obtained for the simultaneous point- 
ing on the first and second observations, respectively. The 
dead time correction fractions for the first ASCA data set are 
76.1% for GIS2 and 77.9% for GIS3, determined by reference 
to the count-rate monitor data with little dead time, using 
the method described in Makishima et al (1996). Similarly, 
for the second observation, these are calculated as 85.6% 
and 87.3% for GIS2 and GIS3, respectively. We only use the 
GIS data (0.7-10 keV) since the SIS is strongly affected by 
pileup. 



In order to take into account the calibration uncer- 
tainties, we add 1% systematic errors to the GIS spectra 
(Makishima et al 1996) and 0.5% systematic errors to the 
PCA. We also include 0.5% systematic uncertainty on the 
HEXTE data to be conservative, though this does not sig- 
nificantly change the fit (see Appendix). We allow Nh to 
be free (described using wabs) as the low energy bandpass 
of the ASCA GIS constrains this parameter. The data from 
different instruments are fitted with same model parameters 
except that relative normalization factors are free. Hereafter, 
we use the normalization of the PCA in the fitting results 
so as to be able to directly compare with the PCA analysis 
of Kubota & Makishima (2004) for the disc spectra. Using 
this we find the observed (i.e. absorbed) 0.7- 200 keV fluxes 
of 7.2 X 10~* erg s"^ cm"^ and 8.3 x 10"* erg s"^ cm"^ 
for the first and second observations, respectively, corre- 
sponding to observed luminosities of 2.2 x 10** erg s~^ and 
2.5 X 10** erg s~^ for isotropic emission. Correcting for ab- 
sorption and the energy band gives a bolometric luminosity 
of ~ 4.0 and 4.7 x 10** erg s~^ i.e. about 35 per cent of the 
Eddington limit. 



4 CONTINUOUS DISC-CORONA 
4.1 Fitting the data 

These VHS spectra are dominated by a strong, steep con- 
tinuum rather than the disc, though the disc can still be 
seen. The continuum shape is complex; fits using a single 
power law (together with the disc, and including a smeared 
edge and Gaussian to phenomenologically model the refiec- 
tion components) are completely statistically unacceptable, 
with xl of 2538/195 and 1718/195, for the first and second 
datasets, respectively. This was also shown by Sobczak et al 
(1999), who had to fit the PCA and HEXTE data separately 
in order to get acceptable fits with a power law continuum. 
However, the continuum curvature is complex: a single ther- 
mal Comptonisation model for the tail (with disc, smeared 
edge and Gaussian line) gives a much improved fit over a 
power law, but the fits are not formally acceptable at xt of 
402.9/193 and 261.5/193, respectively. 

The complex form of the tail is shown in more detail 
by Gierliiiski & Done (2003), who used OSSE data together 
with the PCA and HEXTE to extend the bandpass beyond 
500 keV for times around that of the second dataset used 
here. These clearly showed that the curvature seen in the 
tail cannot be matched by Comptonisation from thermal or 
non-thermal electrons (see Fig. 5 of Gierlinski & Done 2003). 
Instead the contirmum shows features of both thermal and 
non-thermal electrons (see also Zdziarski et al 2001 for simi- 
lar spectra from GRS 1915-H05). Here we follow KD04 and 
approximate this by using the thcomp model (Zdziarski, 
Johnson & Magdziarz 1996; now publically available as a 
local model in XSPEC*) to describe the thermal Compton 
scattering, together with a weak power law of F fixed at 2 
to approximate the higher energy (non- thermal) emission. 
We have checked that our results on the disc temperature 
and luminosity are not dependent on this approximation 

* 

http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/newmodels.html 
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Figure 1. The schematic geometry envisaged for each of the models described in the text. To the right of each geometry is a sketch of 
the associated emissivity as a function of radius for the outer disc (red), inner disc (blue), corona (magenta) and total (black) accretion 
flow. 



for the shape of the high energy spectrum by using the 
full hybrid thermal-nonthermal Comptonisation code EQ- 
PAIR (Coppi 1999) on our final model of the Comptonised 
disc-corona described in §6. 

The seed photons in thcomp are set to be from the 
disc (giving a slightly different low energy rollover to that 
produced by the COMPTT model in the XSPEC general re- 
lease which assumes a Wien distribution) . This is equivalent 
to assuming that the Comptonising corona covers the whole 
disc, as shown schematically in Fig. la. 

Following KD04 we use several different descriptions 
of the reflected emission, to try to assess the effect on the 
derived disc parameters of systematic uncertainties in cur- 
rent reflection models. Firstly we use the Gaussian line 
and smeared edge to give a phenomenological description 
of reflection (SEG). Then we replace these components 
with the physical reflection model incorporated in thcomp 
(REF) which includes both full reflected continuum and self- 
consistent iron line emission for arbitrary ionization state 
of the reflected material (Zycki, Done & Smith 1999). Ion- 
ization is important as the lack of photoelectric opacity in 
ionized material means reflection contributes to the low en- 
ergy spectrum so can affect the derived disc temperature 
and luminosity. 

However, the ionised reflection in THCOMP does not in- 
clude Compton up and downscattering of the reflected emis- 
sion in the disc itself. This is important for highly ionized 
and/or hot discs, and leads to strong distortions on the pre- 
dicted line and edge shape in the reflected emission which 
can affect the spectral fitting, either by distorting the in- 
ferred solid angle of reflecting material, or its ionization or 
amount of relativistic smearing (Ross, Fabian & Young 1999; 
Ballantyne, Ross & Fabian 2001; Done & Gierlihksi 2006). 

Fitting with better ionized reflection models which do 
include the Compton scattering on the line and edge is not 



straightforward as these models are tabulated for an illumi- 
nating power law continuum, which is not appropriate for 
these data. Hence we attempt to quantify the effect of these 
systematic model uncertainties in the reflected emission by 
removing the 5-12 keV region (i.e. the iron line and edge fea- 
tures which are sensitive to the disc Comptonisation) from 
the spectral fits. We then refit the spectra with ionization 
of the reflector fixed at 10^ (Xi2), 10^ (Xi3) and W (Xi4) 
in order to assess how much the derived disc luminosity and 
temperature are dependent on this parameter, and so get an 
estimate of the systematic model uncertainties (as in KD04). 

We tabulate results from all free parameters in these 
fits in Table Al in the Appendix. Uncertainties are given 
for Ax^ ~ 2.7, and derived from a full error scan where 
all other parameters are free to vary. We include the de- 
rived disc luminosity, which is also subject to uncertainties 
in distance and inclination, but we stress that this is used 
only to compare with the disc luminosity derived using the 
same distance and inclination. The refiected fraction is also 
dependent on the assumed inclination, but this is not used 
further in the paper. 

Fig. 2a and b shows each i^F^ spectrum, deconvolved 
and corrected for absorption using the best fit SEG model 
(phenomenological reflection), together with residuals. The 
two spectra are qualitatively similar, with the disc emis- 
sion being strongly Comptonised. The only differences are 
that the uncomptonised disc spectrum is slightly easier to 
identify in the first spectrum than in the second, and the 
spectral index of the Comptonised emission is slightly fiat- 
ter. This immediately sets some constraints on the Comp- 
tonising geometry. The corona will scatter only a fraction 
Cf[l — exp(— r)] where C/ is the covering fraction of the 
corona over the disc and r is the optical depth. Thus if ei- 
ther Cf or T were small then the disc spectrum would be 
above the Comptonised tail. Since it is not, the spectra in 
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Figure 2. Absorption corrected intrinsic spectra for eacii of the two observations described in the text. The left and right hand panels 
show the simultaneous ASCA-RXTE data taken before and after the peak, respectively. Both datasets are deconvolved using the best 
fit phenomenological model with SEG reflection, and the ASCA and HEXTE normalisations are set to match the PCA. The residuals 
are from this fit (SEG in Table Al), the dkbbth fits (SEG in Table A2) and the dkbbfth fits (SEG in Table A3), respectively, showing 
that all these different disc-corona models are a good description of the overall spectral shape. 



Fig. 2a and b show that the majority of the disc emission is 
covered by an optically thick corona (see also the discussion 
in KD04). 

Both spectra also show two clear breaks within the 
bandpass of the data, a high energy downturn at ~ 30 keV 
and a low energy downturn at 1.2 keV and 1.8 keV for the 
first and second datasets, respectively. We stress that the 
low energy downturn seen in these data arc not due to ab- 
sorption, since the spectra in these figures are absorption 
corrected, though it is dependent on the assumption that 
the continuum seed photons are from the disc. Our derived 
column of ~ 0.7 ± 0.2 x 10^^ cm~^ is consistent within un- 
certainties with that of 0.8 ± 0.2 x lO'^'^ cm~^ derived from 
Chandra data (Miller et al 2003) and is slightly less than the 
full column through our Galaxy in this direction (0.9 x 10^^ 
cm~^), as expected. 

The existence of the low and high energy cutoffs within 
the observed 0.7-200 keV bandpass mean that these data 
cover the majority of the accretion power output. Thus the 
total luminosity is well defined (modulo uncertainties in the 
distance and inclination which cancel when comparing to the 
disc luminosity derived from high/soft state data), rather 
than being strongly afi'ected by uncertainties associated with 
extrapolating the spectral models outside of the observed 
range. 



4.2 Derived luminosity and temperature of the 
disc emission 

Fig 3 shows the luminosity and temperature for the disc 
dominated high/soft state observations of KD04 to compare 
with those seen from the disc in the two VHS state spec- 
tra used here. Reconstructing the intrinsic disc emission is 
model dependent: a Comptonising corona will distort the ob- 
served disc spectrum. However, the Comptonised spectrum 
itself contains a low energy rollover, directly imprinting the 
energy of the seed photons on the boosted emission. Thus 
with broad bandpass spectra, where the energy range of the 
expected disc temperature is covered by the observations as 
in these ASCA-RXTE data, the intrinsic seed photon tem- 
perature is recoverable. 

The seed photon luminosity is equally observable, de- 
spite the effects of Comptonisation shrouding the disc. 
Compton scattering increases the photon energy by on aver- 
age H^E/E ~ 4r^(fcTc/meC^). Thus the Comptonised lumi- 
nosity is related to the scattered disc luminosity simply by a 
factor ~ (l-\-/S.E/E) ~ 1.3 for the steep Comptonised spec- 
tra observed here. This scattered disc luminosity should sim- 
ply be added to the observed disc luminosity to estimate the 
intrinsic disc luminosity of {Ldi.sc + Lthcomp) /i-S Tims the 
intrinsic disc luminosity cannot be very much smaller than 
the total luminosity, which is robustly estimated by these 
data. This is shown plotted against observed disc temperar- 
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Figure 3. Luminosity versus temperature plots for tlie disc. Filled circles and filled triangles show the data from the disc dominated 
states, while the open circles show the weakly comptonised very high state (KD04). Filled and open squares and open squares show the 
strongly comptonised very high state data from before and afer the peak, respectively. There are five points for each of these strongly 
comptonised spectra, corresponding to the five different reflection descriptions. The solid lines show L oc T'* for rj^ =59, 100 and 200 km. 
Panel a and b show the reconstructed disc emission from simple modelling of the spectra (Table Al), corresponding to a corona covering 
the whole disc (see Fig. la). Panel a plots the observed disc temperature against the inferred disc luminosity estimated by Ldisc + Lthcomp- 
This is a slight overestimate since Comptonisation boosts the seed photon luminosity by a factor ~ 1.3 for these data. Panel b shows a 
more careful reconstruction of the disc luminosity using the photon numbers and inclination of the slab. Either way gives similar results, 
showing that the inferred disc luminosity and temperature is quite robustly estimated by these data, and that they imply the inner edge 
of the disc is much larger than that seen in the disc dominated spectra. Panel c shows the reconstructed inner disc temperature and 
luminosity for an inner disc-corona (Fig. lb). This is not very different to that found from the previous fits. By contrast. Panel d shows 
the results for an inner disc whose energetics are coupled to that of the corona (Fig. Ic). The implied change in inner disc radius is now 
much smaller. 



ture in Fig 3a, and plainly shows a marked discrepancy with 
the constant radius/constant colour temperature correction 
observed from the disc dominated spectra. 

A more careful correction for the effects of Compton 
scattering can be done by noting that this conserves pho- 
ton number so the intrinsic disc luminosity can be simply 
estimated by increasing the observed disc luminosity by a 
factor {aNcomp + Ndi3c)/Ndisc where Ncomp and Nmsc are 
the numbers of photons in the thermally Comptonised and 



disc spectral components, respectively, and a are 2 cos i and 
1 for spherical and slab geometries (see Kubota & Mak- 
ishima 2004; KD04). Panel b shows the reconstructed disc 
luminosity assuming a slab geometry and clearly shows that 
the inferred disc is very different to that seen in the high/soft 
state data. 

Thus, while spectral fitting is almost always non-unique, 
the presence of the low and high energy rollover in the ob- 
served bandpass mean that we can make very robust esti- 
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mates for the intrinsic disc temperature and luminosity, giv- 
ing a marked discrepancy with the constant radius/constant 
colour temperature correction observed from the disc domi- 
nated spectra. Disc dominated spectra at the same luminos- 
ity have a temperature of ~ 1 keV, so for a spectrum dom- 
inated by Comptonisation we expect a low energy rollover 
from the seed photon energy at 4fcTi n ~ 4 keV. This is in di- 
rect conflict with the observation of the low energy rollover 
at 1.2 keV or 1.8 keV which shows the seed photon temper- 
ature is < 0.5 keV or < 0.7 keV. The inferred disc is highly 
luminous but with low temperature, implying that the inner 
radius of the disc is a factor ~ 3 higher in the VHS than in 
the disc dominated spectra. 

Plainly, with these models the disc is not consistent with 
a constant inner radius, constant emissivity, constant colour 
temperature correction between the VHS and the disc dom- 
inated state. At least one of these must be changing. An 
obvious change is for the corona to affect the colour temper- 
ature correction. Irradiation will change the vertical tem- 
perature structure of the disc, but this should lead to an 
increasing temperature for the upper layers, so an increased 
colour temperature correction. Yet the data show a lower 
temperature for the inferred luminosity. Similarly, conduc- 
tion between the corona and disc would also increase the 
disc temperature, as would an increase in the stress at the 
inner disc boundary (Agol & Krolik 2000). This robust dis- 
crepancy between the disc structure in the disc dominated 
and VHS motivates the more sophisticated modelling in the 
rest of the paper. 



5 INNER DISC CORONA, SEPARATE 
ENERGETICS 

Here we build a simple model of an inner disc corona, and fit 
it to the data, incorporating the effects of the geometry as 
well as the radiative constraints from Comptonisation. The 
corona is assumed to exist only at small radii (see Fig lb) 
so it removes only the hottest part of the disc spectrum 
by Comptonising it, while not affecting the spectrum of the 
cooler disc at larger radii. The aim is to test whether this 
more physically realistic model still requires the inferred in- 
crease in disc inner radius seen in the previous section. 

5.1 Models 

We assume that the corona is a slab of constant optical 
depth and temperature on top of the disc, extending from 
Tin to Vc, while the disc itself extends from rin to lOOOrin 
(see Fig. lb). Each radius in the disc is assumed to emit 
blackbody radiation, with the standard temperature distri- 
bution of a multicolour DISKBB i.e. T(r) = Ti^{r /n^)'^^* . 
Where the disc is not covered by the Compton slab then its 
emission is unaffected. Otherwise, it is suppressed by a fac- 
tor exp(— r), where r is determined from the Comptonised 
spectrum using the observed F and fcTo, assuming a slab 
geometry. The blackbody emission from the disc at this ra- 
dius is used as the seed photons for the Comptonising cloud, 
whose spectrum is then normalised so that the total num- 
ber of photons in the Comptonised spectrum is equal to the 
number of photons scattered from the blackbody disc emis- 
sion at that radius i.e. 1 — exp(— r). 



A more careful treatment would have to account for 
the different angle dependences of the seed photons and 
comptonised emission - the number of seed photons which 
emerge unscattered is more accurately exp(— r/ cos i), but 
the rmmber which are scattered into this direction is not 
simply 1 — exp(— t/cosz) as photons scattered into a given 
line of sight come from a range of initial seed photon angles 
(see e.g. Pozdnyakov, Sobol & Sunyaev 1983). There is yet 
more complexity as about half the Compton scattered pho- 
tons will be directed down towards the disc. Those which 
are not reflected will be reprocessed, adding to the intrin- 
sic disc emission (Haardt & Maraschi 1993). We neglect all 
these effects in order to be able to approximately calculate 
the spectrum. 

The shape of the spectrum is determined by the pa- 
rameters of the Compton corona (F and kTg), the size of 
this corona relative to the disc rc/nn, and intrinsic inner 
disc temperature, fcTin. The overall normalisation gives the 
apparent disc radius, rin, for the assumed distance and incli- 
nation, as in the diskbb model. Even though the full diskbb 
spectrum is not seen, the model normalisation is that of the 
uncomptonised disk. The model parameters for the disk com- 
ponent are in effect what the disk would have looked like 
without the Comptonising corona, and so can be directly 
compared with the disk normalisation derived from the disk 
dominant high/soft state. 

Fig. 4a shows the results for this model corona parame- 
ters close to those implied by the VHS data i.e. a corona with 
feTe = 10 keV and t = 2 (which corresponds to a photon in- 
dex of the scattered flux of F = 2.3 as the optical depth for a 
slab geometry is half of that which gives the same spectrum 
for a sphere). The inner disc temperature is set to 1 keV to 
match that seen from the disc dominated spectra at similar 
luminosities (see Fig. 3). The radius of the Comptonising 
region is progressively increased over more and more of the 
disc, from rc/nn ~ 1, 1.33, 1.66, 3.33 and 1000. This changes 
the spectrum from a pure disc blackbody, to a more com- 
plex shape. While this does have the same number of free 
parameters as the previous model, it is somewhat more con- 
strained. The normalisation of the disc and Comptonised 
components is not a completely free parameter, but is set 
by Tc/rin and by the Eissumed slab geometry. In the limit 
where the corona covers all the disc then the spectrum has 
the same shape as that of a thcomp model with disc black- 
body seed photons but with the addition of exp(— r) of the 
initial disc blackbody spectrum which is unscattered by the 
corona. 

The spectrum always shows the imprint of the seed pho- 
ton energy at ~ 1 keV. This is because of the physical re- 
quirements of the model, firstly that the seed photons are 
from the inner disc with standard emissivity, and secondly 
that the Comptonised spectrum has the same number of 
photons in it as were removed from the disc emission. A 
standard disc emissivity disc produces most of its luminos- 
ity in the inner regions, and it is these inner regions which 
are covered by the corona. The Compton energy boost is not 
large for these coronal parameters (F = 2.3, kTe = 10 keV) 
so the Comptonised photons, with their clear seed photon 
rollover, add to the emission at similar energies to those of 
the seed photons. Comptonisation conserves photon num- 
ber, so the normalisation of the spectrum is also similar. 
Thus the optically thick corona does intercept most of the 
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Figure 4. An inner disc corona model for Tc = 1 (red, i.e. equivalent to no corona), 1.33 (green), 1.66 (blue), 3.33 (cyan) and 1000 (black 
i.e. equivalent to a corona which completely covers the disc) xr^^ (see Fig. lb). The left panel shows the results for coronal parameters 
close to those inferred from the data, with F = 2.3 and kT^ = 10 keV. This corresponds to r ~ 2 for a slab corona. The right panel 
shows the effect of increasing the coronal temperature to 100 keV. To keep the same spectral index for higher temperature means that 
the optical depth has to decrease to r ~ 0.4. Thus there is always a clear unscattered disc component (unlike the data), and the intrinsic 
disc temperature of 1 keV is seen. 



inner disc emission, but this Comptonisation does not hide 
it as neither the energy nor number of photons are changed 
substantially. 

Figure 4b shows the same series of spectra for increasing 
coronal radius, but with kTc — 100 keV. Though increasing 
the electron temperature increases the energy boost so mak- 
ing it possible to move the photons away from their seed 
energy, to get the same spectral index of F = 2.3 from this 
higher temperature plasma requires a lower optical depth. 
The corona becomes optically thin so the spectrum now is 
clearly always two component (disc plus Comptonisation), 
and the disc is always also clearly at 1 keV. Although only 
the inner disc photons are now boosted far away from their 
seed photon energy, the low optical depth means that not 
enough of them are Comptonised to hide the hottest disc 
components even for a corona that extends over the whole 
disc. 

A comparison by eye of these models with the observed 
spectra in Fig. 2 shows that the observed low energy rollover 
at ~ 1.8 keV strongly constrains the seed photon energy 
even in an inner disc corona model to be substantially less 
than 1 keV. Thus this model is not going to easily change 
the conclusion of the previous section that the inner disc 
temperature is much lower than is seen at this luminosity 
in the disc dominated high/soft states. These models of an 
inner disc corona above a standard emissivity disc still re- 
quire that the disc inner radius in the VHS is substantially 
larger than in the high/soft disc dominated states. 

5.2 Data fitting 

We quantify the argument above by replacing the separate 
DISKBB and THCOMP components in the spectral fits with 
this new inner disc corona model. We include the thcomp 



reflection model in the disc-corona code, to take into account 
illumination of the outer disc by the Comptonised spectrum. 
Again we use the same range of ionized reflection descrip- 
tions as before in order to quantify the uncertainties which 
arise from our modeling of this component. 

First we fit this model to the simultaneous ASCA- 
RXTE spectra with the coronal radius fixed at r-c/rin — 
1000. This model has one less free parameter than the sep- 
arate DiSKBB-f THCOMP fits shown in § 4 because of the spe- 
cific geometry (slab corona) assumed. The observed spectral 
shape of the Comptonised emission (F and kTc) set r for the 
specific geometry, and this fixes the ratio of normalisations 
of disc and Comptonised emission. The fits are significantly 
worse, showing that our assumed specific continuous slab 
corona geometry is not consistent with the data. However, 
the inferred disc temperatures for the underlying standard 
emissivity diskbb spectrum are ~ 0.48 keV and 0.51 keV for 
the first and second data set, respectively. These are much 
the same temperatures as for the previous fits (though with 
worse due to the assumed geometry), illustrating the ro- 
bustness of reconstruction of the seed photon (intrinsic disc) 
temperature 

We then repeat the fits with the coronal radius as a 
free parameter to correspond to the geometry of Fig. lb 
(again with all five different refiection descriptions). While 
this formally gives the same number of degrees of freedom 
as the separate diskbb-|-thcomp fits, the model is more 
constrained by the geometric and physical assumptions. 
Nonetheless, the fits are indistinguishable in fi'om those 
in § 4, strongly supporting an inner disc-corona geometry for 
these VHS spectra. Table A2 gives the results for the under- 
lying (i.e. intrinsic, before Comptonisation) diskbb temper- 
ature and luminosity. These are plotted in the middle panel 
of Fig. 3b, together with the points from the disc dominated 
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spectra (KD04). Plainly, as shown above, the ASCA-RXTE 
spectra still strongly require a very low temperature disc 
which does not fall on the constant radius/constant colour 
temperature correction line derived for the disc dominated 
spectra. Again the inferred apparent inner disc radius is a 
factor 2 larger than that seen in the high/soft disc domi- 
nated states. 



6 INNER DISC CORONA, COUPLED 
ENERGETICS 

The previous model assumed that the disc emission is unaf- 
fected by the presence of the corona. That is, the inner disc 
corona model assumes that while the disc radiates efficiently 
at all radii, the coronal mass flow only radiates within Tc. 
However, this seems rather artificial. The corona dissipates 
some power, and this must derive ultimately from the ac- 
cretion flow. The disc and coronal flows could be kept sep- 
arate in a continuous corona geometry (e.g. Fig. la) with 
some (approximately constant) fraction of the power and/or 
mass accretion taking place via the corona (Zycki et al 1997, 
Haardt & Maraschi 1993; Svensson & Zdziarski 1994). The 
geometry of an inner disc-corona instead implies that there 
is a single accretion flow at large radii, Mdisc.oo which then 
splits at Tc to form the disc and corona (Fig.. Ic). Again, we 
build a model incorporating these concepts, and fit it to the 
data, to test whether it still requires an increased inner disc 
radius. 

6.1 Models 

We use the specfic disc-corona coupling model of Svensson 
& Zdziarski (1994), which assumes most of the mass accre- 
tion rate is still carried by the disc, but where a fraction 
/ of the accretion power is dissipated in the corona. The 
power dissipated in the optically thick disc is discontinuous 
at Tc, where it drops from Mdisc, oo by a factor (1 — /). This 



implies the disc temperature is also discontinuous, dropping 
to (1 — fy^"^ of that expected from Mdisc,aa (see Svensson & 
Zdziarski (1994) calculate the effect on the disc in the limit 
in which most of the mass accretion rate is still carried by 
the disc, but where a fraction / of the accretion power is dis- 
sipated in the corona. The power dissipated in the optically 
thick disc is discontinuous at rc, where it drops from Mdisc,oo 
by a factor (1 — /). This implies the disc temperature is also 
discontinuous, dropping to (1 — fy^^ of that expected from 
Mdisc, oo- This is shown schematically in Fig. Ic. 

We incorporate these ideas into our inner disc-corona 
code to predict the spectrum from a coupled disc-corona 
system. The code starts with the standard emissivity disk 
to produce an initial guess at the seed photon luminosity and 
temperature and calculates the Comptonised flux as before 
(using r and kTc to determine r, and scattering a fraction 
1 — exp(— r) of the disc photons into the Comptonised emis- 
sion). The total Comptonised fiux from all radii is used to 
calculate a (radially averaged) value of / which is used to re- 
calculate the disk emissivity and temperature, acting as seed 
photons for the Comptonising corona, giving a slightly differ- 
ent value of /. The code iterates once more around this loop 
in order to reach a self-consistent solution. The parameters 
of the model are the same as for the uncoupled inner disc- 
corona system described above. Again, the inner disc tem- 
perature is that of the uncomptonised disk i.e. the diskbb 
which would have been seen without the corona present. For 
the uncoupled disk-corona this is simply the underlying disc 
emission, whereas for the coupled disk-corona the underlying 
disc emission is distorted from that of a diskbb by the effect 
of the corona draining energy from the disk. The 'tempera- 
ture' which parameterizes this coupled disk-corona systems 
is the innermost disc temperature, TJ'"*, which would be seen 
in the limit of / = 0, not the inner disc temperature used in 
making the model which is (1 - f^^^Ti^^ 

Fig. 5a and b shows the spectra obtained from such 
models with the same intrinsic disc (T;',"' — 1 keV) and 
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Figure 6. The model spectra for V = 2.3 and kTa = 10 keV for rc/nn = 3.3. The total spectrum (black) is the sum of the contribution 
from the outer disc (red), uncomptonised inner disc (blue line), and comptonised emission (magenta). The dashed blue line shows the 
intrinsic (before Comptonisation) inner disc emission. The left panel shows the uncoupled disc corona, while the right panel shows the 
coupled system. There is a marked decrease in luminosity in the coupled system due to the assumed emissivity (compare the right hand 
panels of Fig. lb and c), but more importantly there is also a marked decrease in the low energy rollover seen in this spectrum. This 
is mostly due to the decrease in emissivity of the inner disc emission. This lowers the normalisation of the comptonised spectrum, so 
increasing the relative contribution of the low temperature, outer disc. 



coronal parameters as for the separate disc-corona shown in 
Fig. 4a (where V = 2.3, kT^ = 10 keV, implying / ~ 0.5 in 
the corona) and b (where V = 2.3, kTe = 100 keV, implying 
/ ~ 0.3 in the corona). 

Fig. 6 illustrates the differences between the coupled 
and uncoupled coronae by showing the individual spectral 
components for a corona with V = 2.3, fcTe = 10 keV and 
'"c/fin = 3.3. Both spectra have the same Mdisk.oo, and as- 
sume the same inner disc radius (so without the corona both 
would give the same disc temperature of Tin = Tj'^* = 1 keV) 
but there is much less total power in the coupled disc-corona 
system. This is as expected, as there is the additional coronal 
emissivity in the uncoupled disc-corona model. More impor- 
tantly, the low energy rollover in the spectrum which is used 
to indicate the seed photon temperature is much lower in the 
coupled disc-corona. The reason for this is not primarily due 
to the reduction in temperature of the inner accretion flow 
under the corona. This is only a factor (1 — f Y^'^ lower than 
Tj'n* i.e. 0.85 keV, so the low energy rollover in the Comp- 
tonised spectra is not significantly smaller (4A;Tin=3.4 keV 
as opposed to 4 keV; KD04, magenta line). The major rea- 
son for the lower energy rollover in the total temperature is 
the reduction in emissivity of the inner disc in the coupled 
corona system. The dotted blue lines show the intrinsic inner 
disc emission, while the solid blue lines show the emergent, 
uncomptonised disc emission. The intrinsic disc emission is 
a factor (1 — /) weaker in the coupled disc-corona, so the 
Comptonised spectrum is also a factor (1 — /) weaker. A 
standard emissivity diskbb has maximum luminosity from 
the inner disc components, hence the Comptonised emission 
(with its high seed photon rollover temperature) is gener- 
ally dominant. For the coupled disc-corona, the reduction in 
emissivity under the corona means that the luminosity of the 



corona can be less than the luminosity of the outer disc. The 
Comptonised emission in the coupled system still contains a 
fairly high seed photon rollover, at (1 — fY^'^T^^^ , but this 
can be masked by the much lower temperature rollover from 
the outer disc emission at ~ rin(rc/rin)-^'"' ~ 0.4 keV. The 
major effect of the coupled energetics is to decrease the frac- 
tion of energy emitted in the inner disc-corona, so enhancing 
the relative importance of the cooler, outer disc emission. 

6.2 Data Fitting 

We refit the data with this coupled disc-corona model, in 
which the disc emissivity changes from that of a standard 
disc. Again we use all 5 different descriptions of refiection 
to quantify the effect of systematic uncertainties in spectral 
modelling. Details of the fits to the data are given in Table 
A3, and again the models give equally as good a description 
of the data as either the uncoupled inner disc/slab corona 
(§ 5), or the continuous disc/corona system (§ 4). 

The most obvious change from the coupled energetics 
is that the inferred inner disc radius is much smaller than 
in the previous fits. This is due to the same combination of 
two effects as given above (§ 6.1) i.e. firstly and most impor- 
tantly the increased relative strength of the low temperature 
emission from the outer (uncomptonised) disc, and secondly 
the slight reduction in temperature of the inner disc due to 
the energy drained from it to power the corona. The right 
hand panel of Fig. 3 shows the the derived luminosity and 
temperature of the disc without its corona, T/^', compared to 
the high/soft state luminosity-temperature data. The VHS 
points now lie much closer to those from the disc domi- 
nated states than before (left and middle panels of Fig. 3) 
though formally the statistical uncertainties show that the 
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20 % increase in disc radius in the VHS is significant (as- 
suming a constant colour temperature correction). However, 
this is probably completely compatible with a constant disc 
radius given the systematic modelling uncertainties e.g. our 
assumptions about the shape of the corona (slab, with con- 
stant r and Te), the form of disc-corona coupling (all M 
through the disc, as opposed to some fraction of the accre- 
tion also taking place via the corona) and neglect of self- 
consistent reprocessing of the illuminating coronal flux by 
the disc. 



7 SUMMARY OF SPECTRAL FITS 

The VHS spectra seen in some high luminosity Galactic 
black holes are plainly strongly affected by Compton scat- 
tering. The most extreme VHS spectra no longer have a 
clear disc component, separate from the Compton scattered 
flux, but instead show a smooth continuum shape. The clear 
observational inference is that the radii which produce the 
majority of the disc luminosity must be covered by optically 
thick material. 

Assuming the inner disc is covered by an optically thick 
corona obviously makes it difficult to observe the properties 
of the inner disc, though the intrinsic disc luminosity and 
temperature can still be reconstructed from broad bandpass 
spectra. The temperature can be determined from the low 
energy rollover in the Compton scattered spectrum, while 
the disc luminosity is related to the Compton scattered lu- 
minosity simply through the energy boost from Comptoni- 
sation. 

We use the two simultaneous ASCA-RXTE (0.7- 
200 keV) observations of the extremely Comptonised VHS 
spectra from XTE J1550 — 564 to investigate the intrinsic 
disc emission. As shown by KD04, this emission is not con- 
sistent with that seen in the disc dominated state from the 
same object. The temperature is much lower than that seen 
during the disc dominated states of the same luminosity. 

We build a simple model of a slab corona over a disc 
to show that this conclusion holds irrespective of whether 
the corona covers the whole disc, or only the inner portion. 
The data strongly require that either the disc radius has 
increased, or the disc emissivity has decreased (or both). 
We sketch each of these possible VHS geometries in Fig. 7a- 
c, respectively. If the disc does truncate, the accretion flow 
should still extend down to the last stable orbit unless there 
is a powerful outflow (wind/jet) at this point. Thus we show 
a continuous corona inward of the truncated disc as being 
the likely geometry in Figs. 7a and c. 

The one common denominator for all the geometries in 
Figs. 7 is that the disc is strongly affected by the Compton- 
ising corona. The coronal flow is not some separate phase 
of accreting material, e.g. a low angular momentum flow. 
Instead it must couple in some way to the disc. 



8 ADDITIONAL CONSTRAINTS FROM THE 
DATA 

There is additional information in the spectral and vari- 
ability properties of these data which can give further con- 
straints on the geometry. Firstly, the overall shape of the 
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Figure 7. Possible geometries for the VHS, with radii for the disc 
and corona in a) and c) taken from the average of the Xil, Xi2 
and Xi3 model fits to the spectrum after the peak with dkbbth 
(Table A2) and dkbbfth (Table A3), respectively. 



Comptonised spectrum is steep (F > 2), so the average en- 
ergy gain is rather small. This means that the seed photon 
luminosity is similar to the luminosity emitted by the elec- 
trons i.e. to the heating rate of the electrons. This is easy to 
arrange in the corona over the disc since the disc provides 
copious seed photons. However, it is much more diflicult in 
the geometry of Fig. 7a, with the suggested continuation of 
the corona over a large range (factor 2) in radius down to 
the last stable orbit. The majority of the gravitational po- 
tential energy is dissipated in a region where there is no disc 
underneath, so the only seed photons for Compton cooling 
are a small fraction of those from the disc at larger radii. 
The combination of less Compton cooling and more elec- 
tron heating would surely drive the coronal temperature to 
much higher values, leading to a much harder spectrum (as 
in the low/hard state). The smaller disc truncation implied 
by the geometry of Fig. 7c is much less of a problem, as the 
inner coronal flow is much smaller, so is energetically much 
less important. 

The QPO behavior can also give some clues to the ge- 
ometry, though these are somewhat ambiguous as their ori- 
gin is not yet well understood. The strongest QPO is the one 
at low frequency (LF QPO: type C) and this increases in fre- 
quency as the spectra make a transition from the low/hard 
state to one of the soft states (e.g. van der Klis 2000). This 
indicates that the characteristic radius picked out by the 
variability is decreasing with decreasing strength of Comp- 
tonisation in the spectrum (e.g. di Matteo & Psaltis 1999). 
Qualitatively this can be explained if the disc is truncated 
in the low/hard state, and that its inner edge moves inwards 
towards the last stable orbit during the transition. The fact 
that the LF QPO is not at the highest observed frequency 
in these spectra then indicates that the disc is still slightly 
truncated in these data, favouring the geometry of Fig. 7c 
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over that of Fig. 7b where the disc radius is already fixed at 
the last stable orbit. 

Wc quantify the extent of disc truncation using the spe- 
cific model of Stella, Vietri & Morsink (1999), in which the 
LF QPO frequency, vqpo oc r""'^^, is set by Lense-Thirring 
precession of vertical perturbations of the innermost edge of 
the disc. The LF QPO's (Type-C) are at 2.38 Hz and 3.90 Hz 
for the 1st and 2n(i observations respectively, while in the 
weak VHS data, these are found at 5.9- 6.5 Hz as Type-B 
QPOs (Remillard et al. 2002). The ratio of QPO frequen- 
cies between those seen in these strong VHS, and the most 
weakly Comptonised VHS is 2.6 and 1.6, implying a change 
in radius in the relativistic precession of about a factor 1.4 
and 1.2 for the first and second dataset, respectively. This 
assumes that the weakly comptonised VHS has disc close 
to the last stable orbit (KD04). These estimates from the 
QPO's are startlingly close to the radii derived from our 
coupled disc- coronal fits, which imply a disc inner radius of 
~ 1.35— l.lx andl.3— l.lx that inferred from the high/soft 
state data. However, the major problem with such models 
is that the LF QPO is seen in the spectrum of the corona, 
not of the disc (van der Klis 2000; Gilfanov, Revnivtsev & 
Molkov 2003; Sobolewska & Zycki 2006). The corona may be 
driven at this frequency by the disc, or the LF QPO might 
represent a mode of the corona itself, in which case there 
are no constraints on the geometry. 

The high frequency (HF) QPO seems initially more con- 
straining as the rapid timescale implies that it originates 
close to the black hole. This almost certainly arises from a 
resonance, as it appears as two separate features in a 3:2 
ratio (see e.g. the review by Remillard 2005). If this is the 
parametric resonance between vertical and radial epicyclic 
modes of a disc in strong gravity then this requires that 
the HF QPO comes from the innermost stable orbit of the 
disc around a high spin black hole (Abramowicz & Kluz- 
niak 2001). This model would also require that the inner 
disc can be seen directly, in apparent conflict with the basic 
model presented here of an optically thick inner disc corona. 
However, we note that the HF QPO's are weak, and are not 
seen in our datasets (Remillard et al 2002). Instead they 
are detected in VHS spectra from this source where the tail 
is not as strong as seen here (Remillard et al 2002; KD04; 
Remillard 2005). It is possible that the HF QPO comes from 
the small fraction exp(— r) of emission from the disc which 
is able to escape directly (uncomptonised), so is absent in 
our data simply because these lines of sight are rarer when 
the corona is more optically thick and/or covers more of the 
inner disc. However, it is also possible that the HF QPO rep- 
resents a resonance of the corona, in which case there are 
no visibility constraints (or constraints on black hole spin in 
the specific model of Blaes, Arras & Fragile 2006). 

The only observation which seems in direct conflict with 
all models of an optically thick inner disc corona is that of 
the strong, extremely broad iron line claimed for the sec- 
ond set of ASCA data used here (Miller et al 2004) . Again, 
a small component, exp(— r) of the highly skewed line from 
the inner disc could be seen directly, but this is much weaker 
than the feature observed. Either the X-ray tail in the VHS 
is not formed by Comptonisation, or the spectral residuals 
around the iron line region do not indicate a highly relativis- 
tic line from the inner disc. 



9 ALTERNATIVES TO COMPTONISATION 

The conflict with the broad iron line motivates us to review 
the evidence for Comptonisation as the emission mechanism 
for the tail seen in these data. The only feasible alternar 
tive appears to be the jet. X-ray emission from the jet in 
the low/hard state has been proposed by e.g. Markoff et al 
(2001), so here we use this as a guide to what might be ex- 
pected from these models in the VHS. Direct synchrotron 
emission might give rise to a power law, or broken power 
law, or cutoff power law tail (Markoff et al 2001; Wilms et al 
2005; Markoff, Nowak & Wilms 2005). In §3 we showed that 
a power law gave a completely unacceptable fit to these data, 
but a broken power law is as good a description of the con- 
tinuum as our thermal+nonthermal Comptonisation model 
(see Appendix). However, such models run into difficulties 
when extrapolated to lower energies as the steep spectrum 
dominates the bolometric luminosity. Fig. 8a shows this for 
the best fit broken power law and disc continuum model 
for the second dataset. The extrapolated X-ray continuum 
overpredicts the simultaneously observed optical V and B 
points (dcrcddcdncd using A^, = 4.75) of Jain et al (1999) 
by a factor of over 1000. 

Clearly the data strongly require that there is a low 
energy break in the continuum between the soft X-ray and 
optical wavebands. Such a break can be produced from jet 
models at the point where the jet synchrotron emission be- 
comes optically thick to self-absorption (e.g. Markoff, Nowak 
& Wilms 2005; Giannios 2005). However, there are some 
constrains on the energy of this break from variability argu- 
ments. The short timescale r.m.s. variability spectra of these 
VHS PCA data show a monotonic increase in fractional vari- 
ability amplitude with energy form 3 50 kcV, consistent 
with pivoting about a fixed point at ~ 2 keV (Gierlinski 
& Zdziarski 2005). There is no obvious reason for jet models 
to pick out such an energy. This contrasts with models in 
which the X-ray continuum is produced by Compton scat- 
tering of disc photons. This naturally produces a break in 
the Comptonised continuum at energies below the seed pho- 
ton energy (see Fig. 8b), and since this is at 1 — 2 keV there 
is no fine tuning required to explain this energy being picked 
out as the spectral pivot. 

Gierlinski & Zdziarski (2005) show that for the second 
dataset used here, the r.m.s. variability spectrum from the 
PCA is well matched by thermal Compton scattering of the 
disc photons (sec also Rodrigues et al 2004; Zdziarski et 
al 2005; Sobolewska & Zycki 2006 for the same results for 
other VHS spectra). The Comptonisation models actually 
underpredict the optical flux, but this is expected as the 
disc model here only includes direct emission, while repro- 
cessing is likely to be an important contributor to the opti- 
cal emission (e.g. Esin et al 2000; Hynes 2005), and optical 
synchrotron emission from the jet may also be present (see 
Turler et al 2004 for a model of the plateau/state C/class x 
state of GRS 1915-1-105 with a turndown in the optical). 

Given the strength of the disc emission in the VHS then 
it is also rather natural to assume that this is the dominant 
contributor to the seed photons. Again this contrasts with 
the low/hard state, where the disc is nmcli weaker, and the 
jet might be expected to be a more important contributor 
(Markoff et al 2001). However, even for the low/hard state, 
the jet models fail to reproduce the detailed shape of the tail 
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Figure 8. Absorption corrected X-ray and optical spectra of tlie second dataset. The left hand panel shows the best fit broken power 
law model for the X-ray tail. The extrapolation of this down to the optical overpredicts the simultaneously measured V and B data by a 
factor of over 1000. There must be a low energy break in the spectrum. This is naturally produced in Comptonisation models where the 
disc is the origin of the seed photons. The right hand panel shows the best fit thermal Comptonisation model for the X-ray tail. This 
underestimates the optical flux, leaving room for a reprocessed flux component from illumination of the disc and/or some nonthermal 
synchrotron emission from the outer jet. 



(Zdziarski et al 2003), so that the current version of these 
jet models now have the majority of the X-ray emission pro- 
duced by Comptonisation of the disc and cyclo-synchrotron 
photons by thermal electrons at the jet base (MarkofF et al 
2005). Since disc Comptonisation appears to be a component 
of the X-ray continuum even in the low/hard state, this pro- 
cess should be even more dominant with the stronger disc 
emission seen in the VHS. 

Thus both the multiwavelength spectral shape and en- 
ergy dependence of the X-ray variability strongly support 
a disc Comptonisation interpretation of the X-ray tail. Wc 
show in §7 that this is not inconsistent with the presence 
of high (or low) frequency QPO's. A hot (and so presum- 
ably geometrically thick) inner disc corona may actually be 
required to produce the strong jet (McKinney 2005) traced 
by radio emission during this outburst of XTE J1550 — 564 
(Corbel et al 2001). The only obvious inconsistency is with 
the interpretation of the iron line in these data given by 
Miller et al (2004). However, unlike the Comptonisation 
interpretation of the tail, which scorns to be the only vi- 
able mechanism, there docs exist an alternative explanation 
for the line shape. Done & Gierlihski (2005) show that the 
ASCA spectrum in the second VHS observation used here 
can be fit equally well by much less extreme line smearing if 
there is also some ionized absorption present. The underly- 
ing iron line width is then reduced, consistent with originat- 
ing from reflection in the outer disc, which allows the inner 
disc to be mostly covered by the corona. 

Thus there are no data which unambiguously rule out 
an optically thick inner disc corona, while there are many 
observations which strongly support its existence. Such a 
corona will inevitably distort the disc spectrum in the ways 
explored in this paper. 



10 CONCLUSIONS 

We investigate the broadband continuum spectra of two 
VHS observations covering 0.8-200 keV. These have a strong 
X-ray tail, contributing more than half of the total bolomct- 
ric luminosity of the source. The disc in these spectra has an 
apparent temperature which is much lower than that seen 
from the same source at the same luminosity in the disc 
dominated high/soft state. We show that this is not simply 
produced by the effect of an optically thick corona hiding 
the emission from the inner disc, but instead requires that 
the disc changes in some way from that seen in the high/soft 
state, either in area (Fig. 7a), emissivity (Fig. 7b) or both 
(Fig. 7c). 

The most obvious change in these spectra is that there 
is now a corona. Energetically, this must be powered ulti- 
mately by the accretion flow, so there is less power avail- 
able for the disc itself. We build a speciflc example of disc- 
corona coupling, where the corona over the inner disc drains 
power from the disc material under it, changing its emissiv- 
ity (Svcnsson & Zdziarski 1996). We show that this reduces 
the amount of disc truncation required, to a level where it is 
probably not significant given the model uncertainties e.g. if 
some fraction of the accreting mass goes through the corona 
as well as some fraction of the power, or if the geometry of 
the corona is more quasi-spherical than slab-like. However, 
we marginally favour the coupled disc-corona system with a 
slightly truncated disc (Fig. 7c) as this can simultaneously 
explain the LF QPO in the relativistic disc precession model. 
However, since there is as yet no clear consensus on models 
of QPO formation, this cannot give unambiguous informa- 
tion on disc radii. Thus Fig. 7b, where the disc area (though 
not its properties!) is unchanged from the disc dominated 
states is also consistent with the spectra, though this leaves 
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the change in LF QPO as an open question. Conversely, a 
large truncation radius with continuous corona (Fig. 7a) is 
ruled out by the steep Comptonised spectrum seen. 

However, we stress that irrespective of the specific 
model, the strong Comptonisation seen in these very high 
state data clearly requires that the underlying disc proper- 
ties are rather different to those of a disc without a corona. 
Unsurprisingly, the data require that the disc structure 
changes in response to whatever changing dissipation mode 
powers the growth of the coronal structure. 

We also stress that these VHS spectra are very similar 
to those seen from several ULX in terms of having a very 
low apparent disc temperature for their luminosity. This has 
been interpreted as evidence for intermediate mass black 
holes in the ULX, of order IOOOMq (MiUer, Fabian & MiUer 
2004). Obviously this is not an option in XTE J1550 - 564, 
where the black hole mass is constrained by dynamical stud- 
ies to be a normal stellar remnant black hole mass of order 
^ IOM0 (Orosz et al 2002). We suggest that the same pro- 
cess, of a corona draining energy from the inner disc, may be 
operating in these ULX. Many of these would then be com- 
patible with fairly normal binary systems, where the black 
hole mass is at the upper end of that which can reasonably 
form from models of stellar collapse (typically 30Mq or up to 
~ 80 Mq from recent calculations of binary evolution; Bel- 
czynski et al. 2004) accreting at high mass accretion rates 
(0.1-few times the Eddington limit). We will explore this 
possibility in depth in a subsequent paper. 
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APPENDIX A: DETAILS OF SPECTRAL FITS 

We use the second dataset to quantify some aspects of the 
fits, with a baseline model of a broken power law tail (to- 
gether with disc and SEG reflection gives xt = 146.3/193 
for kTdiac ~ 0.671q;q5 keV. By contrast, using the newest 
response (pcarsp 10.1) gives xl = 183.3/193 for kTdisc = 
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0.77 ± 0.06 keV. The shift is within the uncertainties, but 
more important is systematic: representative high/soft state 
spectra also have similarly shifted temperatures so that 
the difference between the disc properties in the VHS and 
high/soft state is robust. Replacing the HEXTE data with 
a file without the 0.5 per cent systematic error gives xt = 
147.6/193 and kTdisc = QMt^M keV, nearly identical as 
the statistical uncertainties are all larger than 1 per cent. 
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Table Al. Best fit pcirameters obtained by a simultaneous ASCA/i?XTi? observation. The continuum is always described by a multi- 
colour disc spectrum, plus a low temperature thermal comptonised component, together with a power law extending to high energies. 
The five different fits show the results of the different descriptions of refiection. These are a narrow gaussian line and smeared edge in 
SEG, while the other models include the full reflection spectrum (line and continuum) in the thcomp model. In the case of Xil, Xi2 

and Xi3, the data in 5-12 kcV is excluded from the spectral fit, is fixed to 30 Rg, and tfic value of ^ is fixed to 10^, 10^, 10^, for 
models Xil, Xi2, and Xi3, respectively. Normalization factors with the PCA data are used to calculate luminosities. 



model 


ATh 
102lcm-i 


Tin 

keV 


Tthc 


keV 


^disk 
(A^disk) 


Lthc(slab) 

(iVthc) 


^pow 


smedge,line (a) 
refiection (b-e) 


xVdof 


Before peak 


SEG 


6.5 ±0.2 


45+0 02 


2.08iro^ 


10.0 ±0.5 


1.03 
(37.6) 


3.14 
(58.0) 


0.86 


E = 7.9t0-3 keV 

_ _ n O+0.2 
'max — ^•'-*_o 1 

width= 3+1 keV 
E = 6.4810-03 keV 
EW=81 ±10 eV 


201.4/192 


REF 


5.9 ±0.1 


51+001 
"■^'--0.02 


2.10 ±0.01 


11 6+0-5 


0.92 
(30.1) 


3.05/2.64 
(51.0/44.9) 


0.76 


n/2-n- = 0.23+0;04 
^ = 6.810-8 X 103 

Rin = 94lgg® Rg 


230.2/194 


Xil 


6 5+°-2 
"■^-0.3 


0.45 ± 0.03 


9 97+0.02 
^•^'-0.03 


14 ±1 


0.86 
(31.0) 


3.52/2.95 
(80.2/64.4) 


0.76 


Q/2n = 0.7610-08 
(« = 102) 


146.3/146 


Xi2 




0.50 ±n. 02 


2.14 ±0.01 


11.5 ± 6 


0.90 

(29.8) 


3.14/2.69 

157.0/ IS. 5) 


0.81 


n/27r = 0.37 ±0.07 

U = 10-'' ) 


161.2/146 


Xi3 


6.2 ±0.2 


0.48 ± 0.02 


2.10 ±0.01 


11.3 ±0.6 


0.97 
(33.5) 


3.10/2.41 
(53.5/42.7) 


0.81 


0/27r = 0.4210-08 
= 10*) 


157.8/146 


After peak 


SEG 


6.7 ±0.2 


0.55 ±0.02 


9 ,9 + 0.02 
^••J^_0.03 


12 ± 1 


1.19 
(36.1) 


4.09 
(80.9) 


0.69 


E = 7.5 ± 0.2 keV 
Tmax = 0.6 ± 0.2 
width= 5I2 keV 

E = 6.4 ±0.1 keV 
EW=4ll^2 eV 


157.8/192 


REF 


7.1 ±0.1 


44+0 06 


2.74 ± 0.03 


> 120 " 




(0) 


5.87/5.32 
(163.9/161.9) 


0.29 


n/27r = 1.7+0-^ 
? = 5l? 

Rin = 12 dl 2 Rg 


226.2/194 


Xil 


7 1+0.3 
'■^-0.2 


0.53 ± 0.03 


2.64 ± 0.03 




0.89 
(27.6) 


5.20/4.00 
(139.5/98.5) 


0.46 


n/2Tv = 1.6+02 
(C = io2) ■ 


113.5/146 


Xi2 


6.5 ±0.3 


0.60 ±0.02 


2.42 ± 0.02 


16 ±2 


1.17 
(32.9) 


4.10/3.29 
(80.9/64.4) 


0.62 


0/27r = 0.73lO;0^ 
= 103) 


131.6/146 


Xi3 


"■"-0.1 


55+0 02 


9 QC-+0.02 




1.13 
(34.2) 


4.18/3.16 
(80.8/63.1) 


0.64 


Q/2^ = 0.6210-04 

a = 104) 


127.9/146 



" The best fit value appeared in > 200 keV. 
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Table A2. Same as Table Al but with dkbbth instead of diskbb plus thcomp to describe a simple (uncoupled energetics) inner 
disc-corona. 



model 


lO^lcm-l 


Tin 

keV 


[km] 


Rc/rin 


Tthc 


keV 


disk-corona"^ 


smedge 
reflection 


xVdof 


I'x-Iiirc- peak 


SEG 


6 6+°-^ 
"■°-0.3 


60+° °* 

U.DU_Q 02 


12911° 


1 r;Q+0.03 
j..oa_Q 02 


2.08l°;°i 


9 9+°-* 


L = 

disc 
T = 


4.20 
= 2.79 

: 2.3 


E = S.otoA 

Tmax = 0.29lg;gg 

width = 2.6_j g 
E = 6.48 ± 0.07 
EW = 80 ± 10 [cV] 


204.9/192 


rlJir 


D.i it U.Z 


U.04 it U.Uz 


1 nv+S 
107_g 


1 c;a+0-01 
J^-3o_0.02 


o 1 1 -i_ n no 


ii.<5 It U.4 


T 

disc 
T = 


A r\A 
4.U4 

= 2.48 

: 2.0 


O /0-77- n OO+0.02 

€ = (3.411;;;) X 10-^ 

Rin = isoI^JJq 


OOT Q /I nA 

zz t .0/ iy4 


Xil 


6.7 ±0.2 




1701^5° 


1 -9+0.05 
^•'^-0.04 


2.27 ± 0.03 


14 ±1 


L = 

disc 
T = 


4.32 
= 3.81 
: 1.6 


n/27r = 1.0 ±0.2 


148.9/146 


Xi2 


6.21°;? 




1101?° 


1.57 ±0.03 


2.12 ±0.01 


11 4+0.5 


L = 

^int 
disc 

r = 


4.04 
= 2.46 
: 2.0 


n/27r = 0.48lO;io 


159.6/146 


Xi3 




62+0 02 


iiolf 


1.46 ± 0.03 


2 094+° °°'' 


11 O+0.6 


L = 

^int 
disc 
T = 


3.98 
= 2.31 
: 2.1 


n/27r = 0.7±0.1 


157.0/146 


After peak 


SEG 


6.9 ±0.2 


0.71 ± 0.03 


110l«2 


1.70 ±0.05 


Q QQ+0.02 

^■•5^_0.03 


12 ± 1 


L = 

disk 
r = 


5.37 
= 3.97 
: 1.7 


E = 7.6 ± 0.2 

_ _ t;+0.3 

width = 4+i 
E = 6.39 ± 0.08 
EW = 50 ± 10[eV] 


151.1/192 


REF 


6.6 ±0.1 


72+° °2 


105111 


1 9+°-2 




20 ±1 


L = 

disc 
T = 


5.33 
= 3.83 
: 1.1 


n/27r = 0.64l°-J3 

^ = (1.4l°, ?J X 102 

Rin = 3llg 


228.5/194 


Xil 


6.9 ±0.1 


64+0-04 


130 ± 15 


9 9 + 0.3 
^•^-0.2 


9 c.o+0.02 


20li 


L = 

disc 

r = 


5.48 
= 3.66 

: 1.0 


Q/27r = 1.2l°-2 


127.0/146 


Xil2 




70+0 03 


iiollS 


1.7 ±0.1 


9 40+0.02 


16l? 


L = 

disc 
T = 


5.35 
= 3.75 

: 1.4 


f2/27r = 0.6±0.1 


121.2/146 


XilS 


6.8 ±0.2 


0.66 ±0.03 


120 ± 10 


, co+0.03 

^•a^_o.o6 


9 qckO.008 


16 ± 1 


L = 

disc 

r = 


5.37 
= 3.47 
: 1.4 


n/2TT = 0.9to.l 


124.2/146 



" The unabsorbed bolometric luminosity of the dkbbth, L, is calculated for inclination 70° by integrating the model spectrum by 
0.001-300 keV. The intrinsic disc luminosity, i'dJgj;, is calculated by referring to observed Hn and Tj^ as Anrf^aT^. 
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Table A3. Same as TaMc Al but with dkbbftii instead of dihkdb ])Ius tiicomp to doscribo an irmcr thsc-corona witii conplcd energetics. 



model 


AfH 
lO^^cm-l 


^int 
in 

keV 


[km] 


i?c/n„ 


Tthc 
keV 




disc-corona" 


smedgc 
reflection 


xVdof 


Before peak 


SEG 


6.4 ±0.1 


U.»»_o.04 


7,+i0 




2.07 ±0.02 


9.8 ±0.5 


L = 4.03 
/ = 0.58 
T = 2.3 


E = 8.0+°-^ 

_ _ r,+0.3 
'max — '-'•'-'_Q ^ 

width = 3l? 
ii = 6.48 ± 0.07 
EW = 85 ± 10 


204.7/192 


REF 


6.1 ±0.1 


0.91 ±0.02 




9 n+0.07 

^■•^-0.05 


9 1,1+0.007 
^■-^^^-0.005 


11 6+°-* 
^-^■"-0.4 


L = 3.94 
J — U.5b 
T = 2.0 


n/27T = 0.3110;°'* 

? - (,2.7_o gj X 10 
= 160l>f°° 


230.9/194 


Xil 


6.4 ±0.1 


0.801H! 




2 43+0.06 


9 90+0.02 
^■■'"-0.03 


13.71?:? 


L = 3.73 
/ = 0.51 
T = 1.7 


n/27r = 0.9 ±0.02 


150.3/146 


Xi2 


6.1 ±0.1 


0-891H? 




9 O+0.2 
^■•^-0.1 


9 1 qy + 0.007 


n 5+°* 

11-3-0.8 


L = 3.64 

/ = 0.55 
r = 2.0 


n/27r = 0.49l°;?^ 


159.3/146 


Xi3 


6.1 ±0.1 


0.83in^ 




9 1 , +0.06 
^■^-^-0.10 


9 1 ,-,1+0.002 

Z.lUl_g (,j^4 


11 6+"-^ 


L = 3.43 
/ = 0.56 
T = 2.1 


n/27r = 0.70l°;i!j 


157.9/146 


After peak 


SEG 


6.7 ±0.2 


0.95iH^2 




9 O9+0.06 
^■•^^-0.07 


2.31 ±0.02 


12 ± 1 


L = 5.23 
/ = 0.47 
T = 1.7 


E = 7.6 ±0.2 

'T'max — 0-5^Q 2 

width = 4l^' 
E = 6.39 ± 0.08 
EW = 50 ± 10 


154.8/192 


REF 


6.81°;? 




86^11 




2 48+0-03 
^■''0-0.02 


21.8l?:^ 


L = 5.48 
/ = 0.39 

T = l,i 


n/2TT = 0.7 ±0.1 

f = (1.1 ±0.2) X 102 
"ill 44 


207.2/194 


Xil 


7 9+0.2 


74+0-03 


120 ± 20 


2.9 ±0.3 


9 C7+0.03 
^■^'-0.05 




L = 5.30 
/ = 0.36 
T = 0.97 


Q/27T = l.lfol 


114.8/146 


Xi2 


6 7+°-2 

"■'-0.1 


89+°-°^ 






2.41 ±0.01 


16 ± 1 


L = 4.93 
/ = 0.43 
T = 1.4 


n/27r = 0.62l°;J0 


123.2/146 


Xi3 


6 7+°-2 

"■'-0.1 


0.83 ± 0.03 




1 97+0.06 
-^■o' -0.10 


9 q^+0.03 


17 ± 1 


L = 4.39 
/ = 0.45 
r = 1.4 


n/27r = 0.9lg;? 


128.1/146 



" The unabsorbed bolometric luminosity of the dkbbfth, L, is calculated for inclination 70° by integrating the model spectrum by 
0.001-300 keV. In this fit, the intrinsic disc luminosity calculated as 47rr?j(Tl]'^* 
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